PRDTEINS 


In  silico  analyses  of  substrate  interactions  with 
human  serum  paraoxonase  1 

Xin  Hu,1  Xiaohui  Jiang,1  David  E.  Lenz,2  Douglas  M.  Cerasoli,2  and  Anders  Wallqvist1* 

1  Biotechnology  HPC  Software  Applications  Institute,  Telemedicine  and  Advanced  Technology  Research  Center, 

US  Army  Medical  Research  and  Materiel  Command,  Fort  Detrick,  Maryland  21702 

^  Physiology  and  Immunology  Branch,  Research  Division,  US  Army  Medical  Research  Institute  of  Chemical  Defense, 
Aberdeen  Proving  Ground,  Maryland  21010 


ABSTRACT 


INTRODUCTION 


Human  paraoxonase  (HuPONl)  is  a  serum  enzyme 
that  exhibits  a  broad  spectrum  of  hydrolytic  activ¬ 
ities,  including  the  hydrolysis  of  various  organo- 
phosphates,  esters,  and  recently  identified  lactone 
substrates.  Despite  intensive  site-directed  mutagene¬ 
sis  and  other  biological  studies,  the  structural  basis 
for  the  specificity  of  substrate  interactions  of 
HuPONl  remains  elusive.  In  this  study,  we  apply 
homology  modeling,  docking,  and  molecular 
dynamic  (MD)  simulations  to  probe  the  binding 
interactions  of  HuPONl  with  representative  sub¬ 
strates.  The  results  suggest  that  the  active  site  of 
HuPONl  is  characterized  by  two  distinct  binding 
regions:  the  hydrophobic  binding  site  for  arylesters/ 
lactones,  and  the  paraoxon  binding  site  for  phos- 
photriesters.  The  unique  binding  modes  proposed 
for  each  type  of  substrate  reveal  a  number  of  key 
residues  governing  substrate  specificity.  The  poly¬ 
morphic  residue  R/Q192  interacts  with  the  leaving 
group  of  paraoxon,  suggesting  it  plays  an  important 
role  in  the  proper  positioning  of  this  substrate  in 
the  active  site.  MD  simulations  of  the  optimal  bind¬ 
ing  complexes  show  that  residue  Y71  undergoes  an 
“open-closed”  conformational  change  upon  ligand 
binding,  and  forms  strong  interactions  with  sub¬ 
strates.  Further  binding  free  energy  calculations  and 
residual  decomposition  give  a  more  refined  molecu¬ 
lar  view  of  the  energetics  and  origin  of  HuPONl/ 
substrate  interactions.  These  studies  provide  a  theo¬ 
retical  model  of  substrate  binding  and  specificity 
associated  with  wild  type  and  mutant  forms  of 
HuPONl,  which  can  be  applied  in  the  rational 
design  of  HuPONl  variants  as  bioscavengers  with 
enhanced  catalytic  activity. 
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Serum  paraoxonases  (PONs)  constitute  a  family  of  mammalian 
enzymes  with  three  members  including  PON1,  PON2,  and 
PON3.1  PON1  has  been  intensively  studied  in  past  decades  owing 
to  its  ability  to  inactivate  various  organophosphorus  compounds 
such  as  the  nerve  agent  sarin  and  the  pesticide  paraoxon,  from 
which  its  name  was  derived  historically.  Human  paraoxonase 
(HuPONl)  is  a  355-amino  acid  (43-kDa)  glycosylated  protein 
expressed  mainly  in  the  liver  and  is  exclusively  associated  with 
the  “good  cholesterol”  high-density  lipoprotein  (HDL).^ 
Although  the  physiological  substrate  and  biological  function  of 
HuPONl  are  still  unclear,  studies  have  shown  that  it  is  involved 
in  a  wide  range  of  physiologically  important  hydrolytic  activities 
such  as  drug  metabolism  and  antiatherosclerosis.1 ^  HuPONl 
exhibits  a  broad  spectrum  of  catalytic  activity  in  the  hydrolysis  of 
various  substrates,  including  lactones,  thiolactones,  carbonates, 
esters,  and  phosphotriesters.  Recent  studies  suggested  that 
HuPONl  is  most  likely  a  lactonase,  whereas  the  phosphotriester¬ 
ase  and  esterase  activities  can  be  categorized  as  “accidental”  or 
“promiscuous.”^-^  A  crucial  question  arises  as  to  how  these 
diverse  activities  are  mediated  in  HuPONl  and  what  structural 
basis  and  catalytic  mechanisms  are  associated  with  the  hydrolysis 
of  these  different  substrates. 

The  three-dimensional  structure  of  a  hybrid  mammalian 
recombinant  PON1  variant  obtained  by  directed  evolution 
(rePONl)  has  recently  been  determined,  providing  the  first  struc¬ 
tural  information  on  this  hydrolase  family.  ^  The  structure 
revealed  a  six-bladed  (3 -propeller  scaffold  resembling  the  Loligo 
vulgaris  diisopropyl  fluorophophatase  (DFPase),^  but  with 
unique  additions  of  helices  at  the  top  of  the  propeller  forming  a 
closed  active  site.  Two  calcium  ions  are  located  in  the  central 
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tunnel  of  the  propeller;  one  is  buried  (termed  the  “struc¬ 
tural  calcium”),  whereas  the  solvent- exposed  “catalytic 
calcium”  is  located  at  the  bottom  of  the  active  site.  On 
the  basis  of  substrate  titration  profiles,  a  His-His  dyad 
(HI  15  and  H134)  in  the  active  site  was  characterized  and 
postulated  to  act  as  a  general  base  and  proton  shuttle, 
directly  involved  in  the  catalytic  mechanism.  1^  However, 
other  studies  also  showed  that  the  H115W  mutant  of 
HuPONI  retains  activity  with  paraoxon,  arguing  that 
HI  15  is  important  for  substrate  binding  and  specificity, 
but  does  not  directly  participate  in  catalysis.^  Blum 
et  al .14  studied  the  substrate  binding  to  DFPase  and 
proposed  an  alternative  catalytic  mechanism  where  the 
calcium  coordinating  residue  D229  in  DFPase  acts  as  a 
nucleophile  in  the  catalytic  process.  Since  an  aspartate 
residue  (D269)  in  HuPONI  is  similarly  positioned  for  a 
nucleophile  attack  in  the  PON1  active  site,  this  mecha¬ 
nism  may  also  be  part  of  HuPONl’s  functionality. 

The  substrate  recognition  and  catalytic  reactions  of 
HuPONI  are  considered  to  be  a  complex  process  likely 
mediated  via  multiple  mechanisms  and  residues.  The 
three  activities  of  HuPONI  (phosphotriesterase,  esterase, 
and  lactonase)  have  been  extensively  investigated  through 
kinetic  and  site-directed  mutagenesis  studies  to  delineate 
substrate  specificity.  A  number  of  residues  have  been 
identified  as  important  for  esterase  and  PON  activities, 
in  that  substitution  of  these  residues  results  in  diminu¬ 
tion  or  loss  of  function.  These  residues  include  L69, 
HI  15,  H134,  D169,  F222,  D269,  H285,  F292,  T332, 
V346,  W281,  etc  A  15,1 6  The  PON  activity  of 

HuPONI,  however,  is  much  weaker  than  the  esterase  ac¬ 
tivity.  Further  studies  showed  that  residues  affecting  the 
lactonase/esterase  and  phosphotriesterase  activities  seem 
to  be  located  in  different  regions  of  the  active  site.  1 
Although  the  participation  of  particular  amino  acid  resi¬ 
dues  in  the  activity  of  HuPONI  can  be  inferred  through 
directed  evolution  and  mutagenesis  experimental  studies, 
a  model  of  the  detailed  architecture  of  the  catalytic  bind¬ 
ing  site  at  the  atomic  level  is  essential  for  a  better  under¬ 
standing  of  substrate  specificity  and  catalytic  activity.  For 
instance,  the  HuPONI  R192-polymorphism  yields  a  10- 
fold  increase  in  paraoxon  hydrolysis,  whereas  the  Q192 
form  hydrolyzes  other  substrates,  for  example,  soman 
and  sarin,  more  rapidly  than  the  R 192 -polymorph  wild 
type. I**  The  lactonase  activity  is  largely  retained  in 
HuPONI  variants,  but  phosphotriesterase  and  esterase 
activities  are  sensitive  to  mutations.^  Because  HuPONI  is 
being  considered  for  development  as  a  catalytic  biosca¬ 
venger  of  organophosphorus  nerve  agents,  structural 
insights  into  the  binding  interactions  of  HuPONI  with  sub¬ 
strates  are  important  for  the  design  of  efficient  HuPONI 
variants  with  enhanced  phosphotriesterase  activity. 3,20,21 

In  this  article,  we  performed  a  detailed  computational 
analysis  of  HuPONI  interactions  with  various  substrates 
by  means  of  homology  modeling,  docking,  and  molecu¬ 
lar  dynamic  (MD)  simulations.  HuPONI  is  highly  simi¬ 


lar  in  amino  acid  sequence  to  the  recombinant  rePONl, 
allowing  a  reliable  3D  structural  model  of  HuPONI  to 
be  constructed.  The  binding  modes  of  HuPONI  with 
phosphotriesters,  arylesters,  and  lactones  were  probed  by 
docking  representative  substrates  and  their  analogs  into 
an  ensemble  of  protein  conformations.  The  optimal 
binding  conformation  and  energetics  of  the  predicted 
binding  complexes  with  wild-type  enzyme  and  mutants 
were  further  investigated  by  MD  simulations  and  binding 
free  energy  calculations.  To  our  knowledge,  this  is  the 
first  theoretical  study  at  the  atomic  level  of  HuPONI  in 
the  context  of  various  substrates.  The  results  are  in  good 
agreement  with  experimental  observations,  providing  a 
detailed  molecular  model  of  the  structural  and  energetic 
determinants  of  HuPONI  substrate  binding  and  specificity. 

METHODS 

Substrate  compounds 

Three  types  of  substrates  of  HuPONI  (phospho¬ 
triesters,  arylesters,  and  lactones)  are  studied  in  the  pres¬ 
ent  work.  Table  I  shows  the  representative  substrates 
including  phenylacetate,  4-acetoxybenzaldehyde,  y-butyro- 
lactone,  y-undecanoiclactone,  8-valerolactone,  and  para¬ 
oxon.  The  kinetic  parameters  of  KM ,  /ccat,  and  kcat/KM 
were  taken  from  the  work  of  Khersonsky  and  Tawfik.^ 
These  experimental  data  were  obtained  from  the  rePONl 
variant  G2E6.H  Because  the  enzymatic  properties  and  bi¬ 
ological  activities  of  rePONl -G2E6  are  very  similar  to 
the  wild- type  HuPONI,  and  sequence  variations  between 
rePONl -G2E6  and  HuPONI  are  in  regions  that  do  not 
affect  their  active  sites  and  overall  structured  the  experi¬ 
mental  data  obtained  from  rePONl  are  in  general  valid 
for  the  wild-type  HuPONI.  Three-dimensional  coordi¬ 
nates  of  all  compounds  were  generated  using  Sybyl  7.od 
A  short  minimization  (500  steepest  descent  steps  with 
Tripos  force  field)  was  performed  to  remove  unphysical 
internal  strains  in  the  molecules. 

Homology  modeling 

The  3D  structural  model  of  HuPONI  was  built  using 
the  program  PSPP  package  (Protein  Structure  Prediction 
Pipeline)  developed  in  our  group  (http://www.bioanalysis. 
org/).23  pspp  contains  a  suite  of  open-source  software 
that  predicts  protein  structures  from  sequence  through  the 
integration  of  multiple  programs  including  domain 
boundary  detection,  sequence  homology  search,  fold  rec¬ 
ognition,  homology  modeling,  de  novo  design,  and  model 
evaluation.  For  the  HuPONI  modeling,  the  crystal  struc¬ 
ture  of  rePONl  (PDB  identifier  1V04)  was  used  as  the 
template.^  The  disordered  loop  (72-79)  was  recon¬ 
structed  in  the  modeled  structure  of  HuPON  1  by  applying 
spatial  restraints  and  energy  minimization.  The  resulting 
model  was  optimized  with  an  iterative  approach  until  no 
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Table  I 

Data  Set  of  PON  1  Substrates  Used  in  this  Study 


X.  Hu  et  al. 


Name 


Structure 


Phenylacetate 


4-Acetoxy  benzaldehyde 


y-Butyrolactone 


8-Valerolactone 


y-Undecanoic  lactone 


Paraoxon 


kcaf  (s  ) 
698 


26 


111 


210 


62 


4.8 


Kinetic  parameters  for  substrate  hydrolysis  by  rePON  1  were  taken  from  the  work  of  Khersonsky  and  Tawfik. 


Km  InrM) 
1.2 


2.3 


71 


0.57 


0.60 


0.8 


kcaAMfs-1  M-1) 
595,000 


11,100 


~5300 


370,000 


103,000 


5800 


significant  improvement  was  obtained.  The  3D  fold  of  the 
generated  models  was  verified  by  the  program  PROSA 
II,- ^  and  the  stereochemical  quality  and  protein  structure 
of  the  final  models  were  validated  by  PROCHECK.25 

MD  simulations 

MD  simulations  were  conducted  for  modeled  systems 
in  explicit  solvent  using  the  AMBER  9.0  package  and  the 
Parm99  force  field.26  The  solvated  protein  systems  were 
subjected  to  a  thorough  energy  minimization  before  MD 
simulations  by  first  minimizing  the  water  molecules 
while  holding  the  solute  frozen  (1000  steps  using  the 
steepest  descent  algorithm),  followed  by  5000  steps  of 
conjugate  gradient  minimization  of  the  whole  system  to 
remove  close  contacts  and  to  relax  the  system.  Bond 
lengths  involving  hydrogen  were  constrained  with 
SHAKE^  and  the  time  step  for  all  MD  simulations  was 
set  to  2  fs.  A  nonbonded  cutoff  of  10  A  was  used,  and 
the  nonbonded  pair  list  was  updated  every  25  times 
steps.  Periodic  boundary  conditions  were  applied  to  sim¬ 
ulate  a  continuous  system.  The  particle  mesh  Ewald 
method  was  used  to  calculate  the  long-range  electrostatic 
interactions.-^  The  simulated  system  was  first  subjected 
to  a  gradual  temperature  increase  from  0  to  300  K  over 
100  ps,  and  then  equilibrated  for  500  ps  at  300  K,  fol¬ 
lowed  by  production  runs  of  5-10  ns  length  in  total. 
Constant  temperature  and  pressure  (300  K/l  atm)  were 


maintained  using  the  Berendsen  coupling  algorithm^ 
with  a  time  constant  for  heat-bath  coupling  of  0.2  ps. 
The  resulting  trajectories  were  analyzed  using  the  PTRAJ 
module  from  the  AMBER  package.  The  root-mean- 
square  deviations  (RMSDs)  of  the  backbone  were  calcu¬ 
lated  from  the  trajectories  at  1  ps  interval,  with  the  initial 
structure  as  the  reference.  All  MD  simulations  were  car¬ 
ried  out  at  the  US  Army  Research  Laboratory  Major 
Shared  Resource  Center. 

Docking 

The  AutoDock  4.0  program^  was  applied  for  docking 
the  substrates  to  the  modeled  structure  of  HuPONl. 
Multiple  conformations  associated  with  residue  Y71  and 
the  corresponding  flexible  loop  72-79  were  obtained 
from  MD  simulations  and  used  for  all  docking  studies. 
To  avoid  artificial  steric  hindrance  of  badly  placed  side 
chains  in  the  active  site,  residues  Y71  and  R192  were 
treated  as  flexible  during  the  docking  process  with  Auto- 
Dock  4.0.  For  HuPONl  mutants,  the  corresponding 
mutated  residues  were  also  treated  as  flexible.  The  active 
site  of  the  protein  was  defined  by  a  grid  of  70  X  70  X 
70  points  with  a  grid  spacing  of  0.375  A  centered  at  the 
center  of  mass  of  the  phosphate  ligand  of  the  model 
structure.  The  Lamarckian  Genetic  Algorithm  in  Auto- 
Dock  was  applied  with  50  runs  and  the  maximum  num¬ 
ber  of  energy  evaluations  was  set  to  2  X  106.  Default 
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force  field  parameters  for  calcium  in  AutoDock  4.0  was 
adopted  (r  =  0.99,  8  =  0.55).  The  charge  of  calcium  was 
set  to  +1  to  avoid  overestimating  its  interactions  with 
substrate.  All  the  docking  jobs  were  processed  through 
the  DOVIS  program  developed  in  our  groups  Results 
differing  by  less  than  1.0  A  in  positional  RMSD  of  sub¬ 
strate  were  clustered  and  the  final  binding  conformations 
were  represented  by  the  one  with  the  most  favorable  free 
energy  of  binding.  The  optimal  binding  complexes  were 
subjected  a  stepwise  energy  minimization  and  MD  simu¬ 
lations  in  explicit  solvent  as  described  above. 

Binding  free  energy  calculation 

The  binding  free  energies  were  calculated  using  the 
MM-GBSA  method.32’33  a  set  of  300  snapshots  was 
extracted  from  trajectories  of  binding  complexes  at  10  ps 
intervals  from  the  last  3  ns  of  each  MD  simulations.  The 
molecular  mechanics  interaction  energies  were  calculated 
with  the  SANDER  module  in  the  Amber  9  package.^ 
The  polar  contribution  (Ggb)  was  calculated  using  the 
generalized  Born  model  by  Onufriev  et  ai?^  The  same 
infinite  distance  cutoff  (cut  =  999.0)  was  used  in  all  MM 
and  GB  calculations.  The  nonpolar  contributions  (Gs a) 
were  estimated  using  the  MSMS  algorithm^  according 
to  the  equation:  Gs a  —  y  X  SASA  +  b  kcal  mol-1,  with 
7  and  b  set  to  0.00542  kcal/(mol  A-2)  and  0.92  kcal/mol, 
respectively,  and  the  probe  radius  used  to  calculate  the 
solvent  accessible  surface  area  (SASA)  was  set  to  1.4  A. 
The  entropy  (TS)  was  estimated  through  normal  mode 
analysis  using  the  NMODE  module.  Because  of  the  high 
computational  demand,  30  snapshots  were  used  for  en¬ 
tropy  estimation.  Snapshots  were  selected  at  equally 
spaced  intervals  (100  ps)  of  the  last  3  ns  trajectories  from 
each  simulation  system.  Before  the  normal  mode  analysis, 
thorough  minimizations  were  carried  out  with  a  dis¬ 
tance-dependent  dielectric  constant  (s  =  4r^  where  r^  is 
the  distance  between  atoms  i  and  j)  including  all  non- 
bonded  interactions.  Decomposition  of  the  calculated 
binding  free  energies  was  performed  using  the  Amber  9 
package. 

RESULTS  AND  DISCUSSION 

Three-dimensional  structural 
model  of  HuPONI 

Human  PON1  shares  86%  sequence  identity  to  the 
recombinant  variant  (rePONl).  Given  the  high  homology 
to  the  template  structure,  the  3D  models  can  be  consid¬ 
ered  reliably  constructed.^  Figure  1  shows  a  3D  model 
of  HuPONI  in  complex  with  the  P04-ligand.  Structural 
alignment  of  the  HuPONI  model  to  the  template 
rePONl  exhibited  a  0.23  A  of  RMSD  of  the  backbone 
atoms,  confirming  that  the  fold  is  essentially  the  same. 
The  quality  of  the  modeled  structure  was  also  validated 


Figure  1 

The  modeled  structure  of  HuPONI.  Residues  Y71,  R192,  and  F292,  as 
well  as  the  phosphate  ion  bound  in  the  active  site  are  shown  in  the 
“stick  and  ball”  representation.  The  two  flexible  loops  (LI  and  L2)  are 
shown  in  blue. 


with  PROCHECK.  Ninety-five  percent  of  residues  lie  in 
the  most  favored  regions  of  the  Ramachandran  plot 
(100%  in  allowed  region).  Likewise,  the  Z-scores  com¬ 
puted  by  PROSA  II  are  highly  comparable  to  that  of  the 
crystal  structure  of  rePONl. 

Analysis  of  the  HuPONI  model  revealed  that  structural 
features  observed  in  rePONl,  that  is,  the  six-bladed  13- 
propeller  scaffold,  the  putative  catalytic  dyad  His  115- 
Hisl34  in  the  active  site,  the  calcium-ligating  residues 
centered  in  the  tunnel,  and  the  three  helices  at  the  top  of 
the  propeller  are  well  conserved  in  HuPONI.  The  phos¬ 
phate  ion  modeled  in  HuPONI  was  bound  to  the  cata¬ 
lytic  calcium  in  the  same  manner  as  observed  in  the 
structure  of  rePONl.  Helices  HI  and  H2  exhibit  major 
sequence  variations  between  HuPONI  and  rePONl,  but 
most  of  these  mutated  residues  are  located  on  the  protein 
surface.  The  hydrophobicity  of  the  N-terminal  helix  HI, 
which  is  presumably  associated  with  HDL  binding,  is 
also  retained.  Residues  defining  the  active  site  of 
HuPONI  are  generally  identical  to  those  of  rePONl, 
with  the  exception  of  K192,  which  is  normally  R/Q  in 
HuPONI.  A  previously  reported  model  structure  of 
HuPONI  (PDB  code  1XHR)  was  based  on  a  template  of 
DFPase.^l  In  comparison  with  this  structure,  the  (3-propel- 
ler  folding  motifs  agree  well,  but  the  three  helices  that  are 
essential  for  the  formation  of  the  catalytic  site  are  missing 
from  the  DFPase-derived  model,  resulting  in  significant  dis¬ 
crepancies  among  key  residues  in  the  active  site. 
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Figure  2 

(A)  RMSDs  calculated  for  the  backbone  atoms  of  protein  over  a  time 
course  of  10  ns  for  the  two  simulated  systems.  (B)  The  calculated 
fluctuations  of  CA  backbone  atom  (RMSF)  of  HuPONl  for  the 
simulations  of  the  apo  and  complexed  HuPONl  systems. 


Dynamics  of  HuPONl 

To  gain  insight  into  the  stability  and  dynamic  proper¬ 
ties  of  the  HuPONl  structure,  we  performed  explicit  sol¬ 
vent  MD  simulations  of  both  the  apo  form  and  the 
phosphate  ligand-bound  complex.  Figure  2A  shows  the 
plots  of  RMSDs  calculated  for  the  backbone  atoms  of 
the  two  simulated  systems  over  a  time  course  of  10  ns. 
The  ligand-bound  complex  displayed  a  stable  plateau 
after  2 -ns  of  simulations  with  an  averaged  RMSD  value 
less  than  1.5  A,  whereas  the  RMSD  of  the  apo  form  devi¬ 
ated  more  with  an  abrupt  change  at  ns.  At  this  point 
the  RMSD  increased  from  an  average  of  1.5  to  2.0  A, 
indicating  a  conformational  transition.  To  detect  local 
structure  motions,  the  fluctuations  of  backbone  atom  CA 
were  calculated  from  the  simulation  trajectories.  As 
shown  in  Figure  2B,  both  apo  and  ligand-bound  proteins 
exhibited  significant  fluctuations  within  the  regions 
surrounding  the  active  site,  reflecting  the  flexibility  of 
the  catalytic  binding  site  of  HuPONl.  Except  for  the 


N- terminal  tail  of  HI,  loop  72-79  (designed  as  LI  in  this 
article)  exhibited  the  highest  level  of  fluctuations 
throughout  the  MD  simulations.  This  is  consistent  with 
the  disordered  form  of  this  loop  observed  in  the  crystal 
structure  of  rePONl.  In  addition,  the  H3-extended  loop 
292-300  (designed  as  L2)  also  exhibited  high  fluctuation 
in  the  simulations.  Both  loops  LI  and  L2  are  located  at 
the  front  of  the  active  site,  forming  a  tunnel  entrance  to 
the  binding  pocket.  Therefore,  the  high  flexibility  of  these 
loops  might  play  a  selectivity- determining  role  in  sub¬ 
strate  recognition.4  In  comparison  with  the  apo  form, 
atomic  fluctuations  in  the  ligand-bound  form  were  sig¬ 
nificantly  decreased,  especially  for  loop  LI.  The  motion 
of  a  small  loop  L3  (175-278)  on  the  back  of  the  active 
site  was  found  to  be  highly  correlated  with  loop  L2  in 
the  apo  form,  and  also  stabilized  upon  ligand  binding. 
However,  helix  H2  became  more  mobile  in  the  ligand- 
bound  state,  implying  that  the  dynamic  effects  of  ligand 
binding  are  mediated  differently  across  the  active  site. 

A  detailed  trajectory  analysis  revealed  that  the  active 
site  in  the  apo  form  undergoes  a  “closed-open”  confor¬ 
mational  change  associated  with  the  flexible  loop  LI  and 
residue  Y71.  Figure  3  shows  the  changes  of  the  torsion 
angle  of  the  Y71  side- chain  over  the  time  course  of  the 
simulation,  with  representative  snapshots  extracted  from 
the  trajectories.  The  closed  conformation,  characterized 
by  Y71  pointing  into  the  active  site,  was  predominately 
observed  with  the  phosphate-bound  complex.  It 
remained  stable  over  the  entire  simulation  due  to  the  for¬ 
mation  of  a  hydrogen  bond  between  the  OH  group  of 
Y71  and  the  phosphate  ligand.  In  contrast,  the  closed 
Y71  conformation  in  the  apo  form  was  retained  only  in 
the  early  stage  of  the  simulation.  The  increased  fluctua¬ 
tions  of  loop  LI  in  the  apo  form  caused  a  distinctive 
conformational  change  of  Y71  with  the  side-chain  point¬ 
ing  away  from  the  active  site,  resulting  in  an  open  con¬ 
formation  (Fig.  3B).  The  open  conformation  in  the  apo 
state  was  stabilized  by  the  H -bonding  interaction  between 
the  OH  group  of  Y71  and  the  guanidino  group  of  R192, 
which  seems  to  be  more  energetically  preferred.  A  short 
period  of  transition  from  “open”  to  “closed”  form  was 
also  observed  at  ~9  ns,  indicating  that  the  free  energy 
barrier  between  the  two  conformations  is  rather  low.  The 
flexible  transition  of  Y71,  and  the  plasticity  of  the  active 
site,  may  be  of  importance  for  HuPONl  substrate  bind¬ 
ing,  as  delineated  in  the  sections  below. 


HuPONI/ester  interaction 

With  the  modeled  3D  structure  of  HuPONl,  we 
docked  the  arylester  substrate  phenylacetate  into  the 
active  site  by  taking  into  account  multiple  conformations 
associated  with  loop  LI  and  residue  Y71.  The  most  favor¬ 
able  binding  conformation  predicted  by  AutoDock  is 
shown  in  Figure  4 A.  Several  interesting  features  can  be 
identified  from  this  binding  complex.  The  substrate  binds 
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Figure  3 

(A)  Overlay  of  the  open  and  closed  conformations  of  HuPONI  observed  in  the  MD  simulations.  Residues  Y71,  R192,  and  F292  are  shown  in  sticks 
(yellow  in  the  closed  conformation,  green  in  the  open  conformation).  The  calcium  ion  is  shown  as  a  magenta  colored  sphere.  (B)  The  plots  of  the 
torsion  angles  of  Y71  side-chain  over  the  time  course  of  10-ns  simulations  of  HuPONI. 


in  the  hydrophobic  pocket  of  the  active  site  and  forms 
extensive  hydrophobic  and  aromatic  stacking  interactions 
with  a  number  of  residues  including  1291,  F292,  Y294, 
T332,  and  V346.  Most  of  these  residues  have  been  found 
to  be  associated  with  esterase  activity  of  HuPONI  from 
site-direct  mutagenesis  studies.  ^  The  carboxyl  oxygen  of 
phenylacetate  binds  to  the  catalytic  calcium  atom  in  the 
“oxyanion  hole”  at  a  distance  of  2.81  A.  The  “catalytic 
residue”  HI  15  is  located  on  the  other  side,  pointing  to 
the  substrate  with  a  distance  of  4.38  A  between  the  e-N 
of  imidazole  and  the  carbonyl  carbon  atom.  This  binding 
orientation  seems  to  be  essential  for  the  catalytic  attack 


via  a  water  molecule,  as  proposed  by  Harel  et  a/.,^  by 
mimicking  the  oxyanionic  moiety  of  the  reaction  inter¬ 
mediate  and  stabilizing  the  charged  complex. 

Notably,  residue  Y71  adopts  a  closed  conformation  in 
the  docked  complex,  pointing  into  the  active  site  as  was 
observed  in  the  crystal  structure  of  rePONl.  To  further 
characterize  the  dynamic  behavior  of  phenylacetate  inter¬ 
actions  with  HuPONI,  we  performed  explicit  solvent 
MD  simulations  for  the  bound  complex.  The  results 
showed  that  residue  Y71  underwent  a  conformational 
transition  from  the  closed  form  to  an  open  state;  similar 
to  what  was  observed  in  the  apo  form.  In  Figure  5,  the 


Figure  4 

(A)  The  closed  binding  conformation  of  phenylacetate  with  HuPONI  predicted  by  AutoDock.  (B)  The  intermediate  binding  conformation  of 
phenylacetate  with  HuPONI  observed  in  MD  simulations.  (C)  The  open  binding  conformation  of  phenylacetate  with  HuPONI  observed  in  MD 
simulations.  The  phenylacetate  substrate  is  shown  in  green,  the  side  chain  of  active  site  residues  are  shown  in  sticks  (yellow)  with  backbone  shown 
in  ribbon  (grey),  and  calcium  ion  is  shown  as  a  magenta-colored  sphere. 
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Figure  5 

The  torsion  angel  of  the  side  chain  of  Y71  taken  from  the  5-ns  MD 
simulations  of  HuP  ONI /phenylacetate  binding  complex.  Mode  A 
represents  the  closed  conformation  of  Y71,  and  Mode  C  represents  the 
open  conformation  of  Y71.  Mode  B  indicates  an  “intermediate” 
conformation  between  the  transitions  that  last  for  ~200  ps. 

change  in  the  torsion  angle  of  the  side  chain  of  Y71  indi¬ 
cates  the  opening  of  the  active  site  at  ns.  The  open 
binding  conformation  remained  stable  over  the  rest  of 
simulation  time.  It  is  likely  that  the  open  conformation 
is  more  energetically  preferred  due  to  favorable  binding 
interactions.  As  shown  in  Figure  4C,  phenylacetate  binds 
to  the  hydrophobic  pocket  in  the  same  manner  as  that  of 
the  closed  conformation,  forming  an  aromatic  stacking 
interaction  with  F292.  A  major  difference  is  that  residue 
Y71  adopts  an  open  conformation  compared  with  the 
initial,  closed  state.  As  a  result,  the  open  Y71,  together 
with  L240,  F222,  and  F292,  form  a  hydrophobic  cavity 
for  the  substrate.  In  addition,  the  OH  group  of  Y71 
could  stabilize  the  open  binding  conformation  by  form¬ 
ing  a  hydrogen  bond  with  both  R192  and  the  substrate, 
thus  contributing  significantly  to  the  binding  affinity. 
Trajectory  analysis  showed  that  phenylacetate  bound  in 
the  active  site  also  undergoes  a  significant  conformational 
change  during  the  opening  of  Y71  (see  Fig.  5).  This  “in¬ 
termediate”  binding  conformation,  which  lasts  for  ~200 
ps,  is  characterized  by  the  substrate  phenyl  ring  posi¬ 
tioned  in  a  perpendicular  manner  that  faces  the  side 
chain  of  Y71  and  forms  an  aromatic  stacking  interaction 
(Fig.  4B).  In  comparison  with  the  open  conformation, 
the  “intermediate”  binding  mode  does  not  seem  to  be 
well  accommodated  within  the  hydrophobic  pocket,  as 
the  stacking  and  hydrophobic  interactions  between  the 
substrate  and  residues  F292,  F222,  and  L240  are  mostly 
disrupted. 

Because  phenylacetate  is  positioned  towards  the  en¬ 
trance  of  the  narrow  hydrophobic  pocket  and  thus  pro¬ 
vides  room  for  bulky  substituents,  the  “intermediate” 
conformation  likely  represents  an  optimal  binding  mode 
for  the  analogs  of  phenylacetate  substrates.  To  validate 


this  hypothesis,  we  docked  4-acetoxybenzaldehyde  into 
the  active  site  and  analyzed  the  energetically  favorable 
binding  conformations.  As  expected,  the  binding  of  this 
analog  to  the  hydrophobic  region  is  similar  to  the  “inter¬ 
mediate”  binding  mode  of  phenylacetate  (see  Fig.  6).  The 
phenyl  ring  faces  the  aromatic  side  chain  of  Y71,  allow¬ 
ing  the  acetoxyl  substitute  to  protrude  from  the  hydro- 
phobic  binding  pocket.  The  para-substitute  of  the  ace¬ 
toxyl  group  may  cause  steric  hindrance  within  the 
pocket,  especially  with  the  “lid”  formed  by  residue  F292, 
preventing  binding  more  deeply  in  the  pocket.  MD  simu¬ 
lations  revealed  that  the  4-acetoxybenzaldehyde  binding 
complex  remains  stable  and  no  significant  conformation 
changes  were  observed. 

HuPONI/lactone  interaction 

Recent  studies  strongly  suggest  that  lactones  are  the 
natural  substrate  of  PON1,  with  kinetic  behavior  differ¬ 
ent  from  other  substrate  classes.^  The  rate  of  hydrolysis 
of  lactones  is  substantially  dependent  on  the  KM  values, 
whereas  this  is  not  the  case  for  esters  and  paraoxon.  For 
a  six-member  ring  containing  8-valerolactone,  the  KM  is 
~50  times  lower  than  that  of  a  five-member  ring  con¬ 
taining  y-butyrolactone.  Ring  substituents  of  y-butyro- 
lactone,  such  as  y-undecanoiclactone,  increase  the  bind¬ 
ing  affinities  proportionally  to  the  length  of  the  aliphatic 
substituents,  resulting  in  dramatic  difference  in  kcat/KM 
values.^  To  probe  the  underlying  atomic  nature  of  these 
substrate  binding  interactions,  we  docked  various 
lactones  to  the  modeled  structure  of  HuPONl  and 


Figure  6 

The  binding  conformation  of  4-acetoxybenzaldehyde  with  HuPON  1  as 
predicted  by  AutoDock.  The  substrate  is  shown  in  green,  the  side  chain 
of  active  site  residues  are  shown  in  sticks  (yellow)  with  backbone  shown 
in  ribbon  (grey),  and  the  calcium  ion  is  shown  as  a  magenta-colored 
sphere. 
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Figure  7 

The  binding  conformation  of  7-butyrolactone  (A),  8-valerolactone  (B),  and  y-undecanoiclactone  (C)  with  HuPONI  as  predicted  by  AutoDock.  The 
substrate  is  shown  in  green,  the  side  chain  of  active  site  residues  are  shown  in  sticks  (yellow)  with  backbone  shown  in  ribbon  (grey),  and  the 
calcium  ion  is  shown  as  a  magenta-colored  sphere. 


performed  MD  simulations  of  the  bound  complexes.  The 
binding  modes  of  7-butyrolactone,  8-valerolactone,  and 
y-undecanoiclactone  are  shown  in  Figure  7.  Similar  to 
the  arylester  phenylacetate,  the  small  lactone  substrates 
y-butyrolactone  and  8-valerolactone  are  well  accommo¬ 
dated  in  the  hydrophobic  pocket  of  the  active  site,  imply¬ 
ing  that  lactone  binding  is  also  driven  by  hydrophobic 
interactions.  The  five/six-member  ring  of  y/8-lactones  is 
positioned  in  such  a  way  that  it  is  perpendicular  to  the 
phenyl  ring  of  F292,  forming  stacking  interactions  with 
the  open  form  of  Y71  [Figs.  7(A,B)].  The  long  aliphatic 
substituents  of  y-undecanoiclactone  protrudes  from  the 
pocket  into  the  hydrophobic  tunnel,  participating  in 
extensive  hydrophobic  interactions  with  residues  along 
the  tunnel  wall,  for  example,  Y294  and  F347  (Fig.  7C). 

Unlike  the  phenylacetate  substrate  bound  in  HuPONI, 
the  perpendicular  binding  orientation  of  the  five-  or  six- 
member  ring  would  allow  the  alkoxy-atom  to  point 
upwards  and  form  hydrogen  bonds  with  residues  N224 
and  N168.  MD  simulations  of  the  two  complexes  showed 
that  the  binding  conformation  of  8-valerolactone  is  quite 
stable  with  the  hydrogen  bonding  interactions  main¬ 
tained.  In  contrast,  y-butyrolactone  exhibited  large  fluc¬ 
tuations  with  the  5 -member  ring  flapping  up  and  down 
intermittently  during  the  simulations.  In  fact,  besides  the 
intrinsic  properties  and  the  size  of  the  ring  involved  in 
the  hydrophobic  interactions,  the  higher  experimental 
binding  affinities  of  8-valerolactone  over  y-butyrolactone 
(0.59  vs.  21  mM  for  KM ,  respectively) ^  are  probably  due 
to  the  well-defined  H -bonding  network  with  the  cal¬ 
cium-ligating  residues.  Notably,  with  the  addition  of  ring 
substituents,  the  stability  of  the  binding  interactions  of 
y-undecanoiclactone  was  significantly  increased.  This  is 
evident  from  the  enhanced  binding  affinities  of  y-unde- 
canoiclactone  and  lack  of  steric  hindrance  along  the 


hydrophobic  tunnel  in  the  predicted  binding  complex 
(Fig.  7C). 

HuPON  1/phosphotriester  interaction 

Figure  8A  shows  the  predicted  binding  mode  of  para- 
oxon  in  HuPONI.  Several  distinguishing  binding  features 
can  be  identified  with  the  optimal  paraoxon/HuPONl 
interaction.  First,  paraoxon  is  orientated  in  a  different 
region  in  the  active  site,  opposite  to  the  hydrophobic 
binding  pocket  used  by  ester  and  lactone  substrates.  This 
paraoxon  binding  site  is  mainly  formed  by  helix  H2  and 
loop  LI,  with  the  postulated  catalytic  residues  HI  15  and 
HI 34  located  at  the  bottom.  A  number  of  residues  within 
this  site,  such  as  R192,  S193,  and  F222,  are  known  to  be 
involved  in  PON  activity.  13>21  Second,  the  paraoxon  sub¬ 
strate  is  situated  on  top  of  the  proposed  catalytic  histi¬ 
dine  dyad,  while  the  phosphate  group  is  pointing  to  the 
catalytic  calcium  atom  facing  residue  D269.  Such  a  bind¬ 
ing  mode  is  consistent  with  the  recently  proposed  cata¬ 
lytic  mechanism  in  which  D269  serves  as  the  catalytic 
residue  for  nucleophile  attack  in  the  hydrolysis  of  phos- 
photriester  substrates,  rather  than  the  His  dyad  used  by 
esters  and  lactones.^  Third,  the  leaving  group  of  the 
nitrophenyl  moiety  is  directed  towards  residue  R192, 
potentially  engaged  in  ionic  and  hydrogen  bonding  inter¬ 
actions  with  the  positively  charged  guanidino  group.  This 
interaction  provides  a  plausible  rational  for  the  well- 
established  difference  in  activity  of  the  polymorphisms 
R192Q,  wherein  the  R  has  significantly  increased  PON 
activity  compared  with  the  Q  form.^  By  this  model,  the 
binding  affinity  of  the  nitrophenyl  group  of  paraoxon 
with  the  positively-charged  R192  would  be  much  stron¬ 
ger  compared  to  that  with  neutral  Q192.  The  estimated 
binding  free  energies  of  paraoxon  with  the  two  R192  and 
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Figure  8 

The  binding  conformation  of  paraoxon  with  wild-type  HuPONl  (A),  H115W  mutant  (B),  and  L69V/V346A  mutant  (C)  as  predicted  by  AutoDock. 
The  substrate  is  shown  in  green,  the  side  chain  of  active  site  residues  are  shown  in  sticks  (yellow)  with  backbone  shown  in  ribbon  (green).  Mutated 
residues  (W115,  V69,  and  A346)  are  shown  in  brown.  The  calcium  ion  is  shown  as  sphere  in  magenta. 


Q192  variants  as  calculated  with  AutoDock  are  —3.69 
and  —3.40  kcal/mol,  respectively.  More  significantly,  the 
favorable  binding  interaction  with  R192  would  facilitate 
the  proper  positioning  of  the  substrate  in  the  active  site. 
This  is  critical  to  the  PON  activity,  as  steric  hindrance 
and  inadequate  positioning  of  the  phosphotriester  sub¬ 
strates  seem  to  be  major  factors  for  their  observed  low 
kcat  values.^  Finally,  residue  Y71  adopts  a  more  open 
conformation  in  the  binding  complex  with  paraoxon. 
This  is  essential  for  paraoxon  bound  in  the  active  site  of 
HuPONl,  as  the  closed  form  of  Y71  would  cause  signifi¬ 
cant  steric  interference  with  the  substrate.  MD  simula¬ 
tions  and  binding  free  energy  calculation  indicate  that 
Y71  plays  an  important  role  in  paraoxon  binding,  and 
makes  significant  contribution  to  the  estimated  binding 
free  energies  (discussed  below). 

A  comparison  of  paraoxon  and  phenylacetate  docked 
to  HuPONl  mutants  provides  structural  insights  into  the 
relative  activities  of  HuPONl  with  each  substrate.  Experi¬ 
ments  show  that  the  F222Y  mutation  abolished  the  PON 
activity  of  HuPONl  but  retained  the  esterase  activity, 
with  a  1.5-fold  increase  in  KM  for  phenylacetate.  How¬ 
ever,  replacing  phenylalanine  with  aspartate  (F222D) 
abolished  both  esterase  and  PON  activity.  13>21  As  seen 
in  Figure  8A,  F222  is  located  on  the  top  of  the  binding 
pocket  and  points  to  the  paraoxon  binding  site.  When 
paraoxon  binds  to  the  F222Y  mutant,  the  substituted  res¬ 
idue  apparently  causes  steric  hindrance  and  prevents  the 
substrate  from  inserting  into  the  catalytic  site.  In  con¬ 
trast,  for  the  phenylacetate  substrate,  no  significant  steric 
hindrance  is  predicted.  The  elimination  of  both  esterase 
and  PON  activities  in  the  F222D  mutant  is  consistent 
with  the  mutated  aspartate  residue  creating  an  unfavora¬ 
ble  electrostatic  interaction  with  the  substrates  and  desta¬ 
bilizing  the  positively-charged  binding  center.  Although 
the  effects  of  mutations  at  F222  on  lactonase  activity 


have  not  been  reported,  we  postulate  that  this  activity 
would  also  be  abolished  with  the  F222D  mutant,  whereas 
the  effect  would  be  less  significant  for  the  F222Y  mutant. 

To  further  investigate  the  impact  of  amino  acid  substi¬ 
tutions,  we  modeled  the  binding  interactions  of  paraoxon 
with  HuPONl  mutants  H115W,  L69V,  and  V346A.  The 
PON  activity  of  these  mutants  has  been  extensively  stud¬ 
ied,  and  the  results  showed  a  ~2  fold  increase  for 
H115W,  and  ~4-16  fold  higher  PON  activity  for  the 
V346A  and  L69V  mutants.  H>21  The  predicted  binding 
modes  of  paraoxon  with  H115W  show  that  the  mutated 
residue  adopts  a  conformation  with  the  indole  side  chain 
parallel  to  the  nitrophenyl  group  of  paraoxon,  thus  form¬ 
ing  a  cluster  of  strong  aromatic  stacking  interactions 
(Fig.  8B).  For  the  mutant  L69V/V369A,  the  docked  com¬ 
plex  suggests  that  the  binding  of  the  ethyl  group  of  para¬ 
oxon,  which  points  to  the  bottom  of  the  active  site, 
seems  to  be  energetically  favored  by  the  enlarged  hydro- 
phobic  pocket  formed  by  the  mutated  residues  A346  and 
V69  (Fig.  8C).  The  increased  binding  affinities  of  para¬ 
oxon  to  the  mutant  enzymes  are  also  mimicked  by  the 
calculated  AutoDock  binding  free  energies,  as  being 
—4.54  kcal/mol  for  H115W,  and  —4.19  kcal/mol  for 
L69V/V346A,  when  compared  with  the  wild-type  binding 
energy  of  —3.69  kcal/mol. 

Binding  free  energy  decomposition 

To  gain  further  insight  into  the  structural  basis  and 
energetics  of  substrate  interactions,  we  calculated  binding 
free  energies  using  the  MM-GBSA  method,  and  decom¬ 
posed  energetic  components  into  individual  ligand-resi¬ 
due  pair.32>33  This  allowed  us  to  estimate  the  contribu¬ 
tion  of  each  residue  to  the  binding  affinity,  and  to  pin¬ 
point  residues  important  for  determining  substrate 
specificity.  The  calculated  binding  free  energies  for  the 
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Table  II 

Binding  Free  Energies  of  HuPONI  (R192)  with  Substrates  Phenylacetate,  y-Undecanoiclactone,  and  Paraoxon  Calculated  by  MM-GBSA 


Substrate 

AFelee 

AEvdw 

A^nonpol 

A^Jpolar 

-TAS 

AGpred 

A^exp 

Phenylacetate 

-11.19  (2.81) 

-20.99  (2.16) 

-3.12  (0.11) 

18.80  (2.68) 

13.92  (7.61) 

-2.57 

— 4.39a,  -3.98b 

7-Undecanoiclactone 

-4.25  (2.67) 

-25.00  (3.36) 

-3.97  (0.26) 

12.61  (2.88) 

15.60  (5.86) 

-5.01 

-4.39b 

Paraoxon 

-46.49  (6.78) 

-35.33  (3.57) 

-4.91  (0.12) 

58.58  (4.09) 

18.93  (7.45) 

-9.23 

-5.188,  -4.22b 

The  predicted  free  energies  (AGpred)  are  the  sums  of  the  contributions  (A£elec,  AEvdw,  AGnonpoi,  AGpoiar)  and  —  TAS).  All  values  are  given  in  kcal/mol.  The  standard  devi¬ 
ations  are  shown  in  parenthesis. 

Experimental  data  (KM)  by  HuPONI  (Q192)  taken  from  the  work  of  Yeung  and  CerasoliU 
hKM  by  rePONl  taken  from  the  work  of  Khersonsky  and  Tawfik.^ 


three  substrates  are  shown  in  Table  II,  whereas  Figure  9 
depicts  the  binding  free  energy  contributions  of  key  resi¬ 
dues  identified  in  terms  of  electrostatic,  van  der  Waals 
interactions,  and  changes  in  solvation.  Overall,  the  van 
der  Waals  interactions  and  the  nonpolar  part  of  solvation 
free  energies  are  predominately  favorable  for  the  binding 
of  phenylacetate  and  lactones,  whereas  for  paraoxon  both 
van  der  Waals  and  electrostatic  interactions  contribute 
favorably  to  the  binding  energies.  These  results  are  con¬ 
sistent  with  our  previous  binding  analysis,  in  which  the 
more  favorable  electrostatic  interactions  observed  in  para¬ 
oxon  binding  are  mostly  due  to  the  interactions  with 
R192,  as  revealed  by  the  energy  decomposition  shown  in 
Figure  9.  The  estimated  binding  free  energies  for  the  three 
substrates  are  generally  in  fair  agreement  with  the  experi¬ 
mental  data  that  was  obtained  from  recombinant  Q192- 
HuPONl  and  the  rePONl  variants  (Table  II).  The  pre¬ 
dicted  binding  affinities  are  expected  to  be  higher  than 
those  for  Q192-HuPONl  because  the  variant  used  in  our 
calculation  is  R192-HuPONl,  and  experimental  data  show 
a  ~  10- fold  increase  of  paraoxon  activity  for  R192.^ 

A  number  of  residues  important  for  substrate  interac¬ 
tion  were  identified  from  the  binding  free  energy  decom¬ 
position.  With  a  free  energy  cutoff  of  0.1  kcal/mol  to 
indicate  binding  of  a  residue  to  the  substrate,  18  residues 
are  predicted  to  participate  in  phenylacetate  binding,  16 
in  lactone  binding,  and  24  in  paraoxon  binding  (see  Fig. 
9).  These  “binding  epitopes”  associated  with  the  binding 
of  each  substrate  can  generally  be  categorized  as  the  first 
shell  calcium- ligating  residues  (E53,  D269,  N168,  N224, 
and  N270),  hydrophobic  pocket-forming  residues  (L240, 
V346,  F292,  1291,  F222,  L69,  HI  15,  H285,  T332,  L267, 
and  F347),  and  paraoxon-binding  residues  from  helix  H2 
and  loop  LI  (Y71,  P72,  R192,  L191,  M196,  D188,  and 
PI 89).  Most  of  these  residues  have  been  extensively  stud¬ 
ied  with  mutagenesis  experiments,  showing  a  range  of 
effects  on  substrate  specificity.  For  phenylacetate  and  y- 
undecanoiclactone,  the  majority  of  residues  contributing 
to  binding  energies  lie  in  the  hydrophobic  pocket.  Resi¬ 
due  L240  contributes  —1.24  kcal/mol  to  phenylacetate 
binding,  whereas  F292  contributes  —2.08  kcal/mol  to  y- 
undecanoiclactone  binding.  Both  are  key  residues  posi¬ 
tioned  on  top  of  the  substrate  as  observed  from  the  pre¬ 
dicted  binding  complexes  (Figs.  4A  and  7C).  For  para¬ 


oxon  binding,  Y71  is  the  largest  contributing  residue 
(  —  2.71  kcal/mol),  followed  by  the  calcium-ligating  resi¬ 
due  E53  (  —  2.15  kcal/mol)  and  R192  (  —  2.07  kcal/mol). 
These  residues  form  the  major  interactions  with  the  phe¬ 
nyl  ring,  the  phosphate  group,  and  the  nitro  group  of  the 
paraoxon  substrate  (Fig.  8A).  It  is  worth  noting  that  Y71 
is  also  a  large  contributor  to  the  binding  of  lactones 
(  —  1.95  kcal/mol),  and  for  interactions  with  phenylacetate 
(  —  0.75  kcal/mol).  These  results  are  consistent  with  our 
binding  mode  analysis,  reiterating  the  potential  impor¬ 
tance  of  Y71  in  substrate  interactions.  Interestingly,  the 
first  shell  calcium-ligating  residues  contribute  more  to 
the  binding  of  phenylacetate  and  paraoxon  than  to  lac¬ 
tone  binding. 

As  an  alternative  to  the  experimental/computational  al¬ 
anine  scanning  method, binding  free  energy  decompo¬ 
sition  provides  a  means  to  rapidly  estimate  mutational 
effects  and  guide  the  experimental  design  of  site-directed 
mutagenesis  efforts.  As  discussed  above,  residue  Y71 
seems  to  be  a  candidate  for  experimental  mutation.  Use¬ 
ful  information  may  also  be  gained  from  examining  the 
contributions  of  residues  in  terms  of  electrostatic  interac¬ 
tions,  van  der  Waals  interactions,  and  solvation  energy 
contributions.  For  example,  for  those  residues  predicted 
to  have  unfavorable  electrostatic  interactions  with  para¬ 
oxon,  for  example,  PI 89,  T332,  N270,  D188  (Fig.  9  and 
supporting  Tables  S1-S3),  mutations  to  nonpolar  resi¬ 
dues  could  be  considered.  D188  is  a  surface  residue 
located  on  H2  distal  from  the  catalytic  site.  It  is  poten¬ 
tially  involved  in  paraoxon  binding  through  interaction 
with  the  nitrophenyl  group.  Hence,  mutation  of  this  resi¬ 
due  to  lysine  or  arginine,  together  with  the  neighboring 
PI 89  mutation  (e.g.,  to  alanine,  leucine,  or  arginine), 
could  effectively  improve  paraoxon  binding.  The  binding 
free  energy  can  further  be  analyzed  in  terms  of  energy 
contributions  from  both  backbone  and  side  chain  atoms. 
As  shown  in  Figure  10,  the  side  chain  of  D183  seems  to 
disfavor  paraoxon  binding  due  to  its  charged  nature  and 
possible  electrostatic  repulsion  with  the  phosphate  group 
of  paraoxon.  Therefore,  the  mutation  D183N  might 
improve  the  binding  affinity  of  HuPONI  for  organo- 
phosphate  substrates.  Similarly,  the  side  chain  of  SI 93  is 
not  predicted  to  make  any  contribution  to  the  binding 
energy,  and  likely  generates  unfavorable  polar  solvation 
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Figure  9 

The  substrate  binding  free  energy  contributions  of  key  active- site 
residues  of  HuPONl  for  phenylacetate  (A),  y-undecanoiclactone  (B), 
and  paraoxon  (C).  The  residues  were  selected  based  on  an  interaction 
energy  cutoff  of  0.1  kcal/mol.  The  residue  binding  free  energies  are 
decomposed  into  individual  terms  based  on  the  MM-GBSA  calculations 
(see  supporting  materials  for  further  details).  [Color  figure  can  be 
viewed  in  the  online  issue,  which  is  available  at  www.interscience. 
wiley.com.] 


energies  (Supporting  Table  S3).  This  is  consistent  with 
the  experimental  observation  that  the  S193P  mutant 
increases  phosphotriesterase  activity.  H 

Implications  for  the  catalytic  mechanism 

Figure  1 1 A  depicts  the  overall  architecture  of  the  active 
site  of  HuPONl  in  the  context  of  substrate  binding  inter¬ 
action.  It  is  evident  from  our  modeling  analysis  that  the 
binding  site  is  characterized  by  two  distinct  regions  on 
the  sides  for  substrate  interactions:  the  hydrophobic 
binding  site  for  arylesters/lactones,  and  the  paraoxon 
binding  site  for  phospho  triesters.  A  well-defined  hydro - 
phobic  pocket  is  formed  on  the  side  of  helix  H3,  mainly 
consisting  of  residues  1291,  F292,  Y294  from  the  H3- 
extended  loop  L2,  residues  L240,  F222  on  the  top,  and 
T332,  V346  at  the  bottom  of  the  active  site.  This  hydro- 
phobic  region  is  typically  associated  with  the  lactonase 
and  esterase  activities  of  HuPONl  from  mutagenesis 
studies.  On  the  other  hand,  the  opposite  side  is  pri¬ 
marily  composed  of  polar  residues,  for  example,  R192 
and  SI 93  from  helix  H2,  and  Y71  from  the  flexible  loop 
LI.  Experiments  have  shown  that  PON  activity  is  sensi¬ 
tive  to  amino  acid  substitutions  in  this  region.  17 

Our  predicted  models  for  substrate  binding  support 
the  proposed  catalytic  mechanisms  of  HuPONl  associ¬ 
ated  with  each  type  of  substrate. 17  pigUre  up 
shows  the  superposition  of  the  three  substrates  (phenyl¬ 
acetate,  8-valerolactone,  and  paraoxon)  bound  in  the 
active  site  of  HuPONl.  Except  for  the  difference  in  bind¬ 
ing  orientation  of  the  phenyl  and  six-member  ring  (par¬ 
allel  vs.  perpendicular),  phenylacetate  and  8-valerolactone 
bind  in  the  hydrophobic  site  in  the  same  manner  with 
the  carbonyl  group  pointing  towards  the  proposed 


Figure  10 

The  residue-based  binding  free  energy  contributions  of  the  backbone 
(red)  and  side  chain  (blue)  for  paraoxon  bound  HuPONl  (see 
supporting  Tables  S1-S3  for  further  details).  [Color  figure  can  be 
viewed  in  the  online  issue,  which  is  available  at  www.interscience. 
wiley.com.] 
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Figure  11 

(A)  A  schematic  representation  of  the  overall  architecture  of  the  active 
site  of  HuPONI.  The  ester/lactone  binding  region  (in  cyan)  and 
paraoxon  binding  region  (in  yellow)  are  shown  as  circles,  and  the 
calcium-ligating  center  with  postulated  catalytic  residues  His  dyad  and 
D269  are  shown  in  the  red  region.  Active-site  associated  helices  (H2 
and  H3)  and  loops  (LI  and  L2),  as  well  as  residues  (R192,  F292,  and 
Y71),  are  also  illustrated.  (B)  Superposition  of  three  substrates 
phenylacetate  (cyan),  8-valerolactone  (green),  and  paraoxon  (yellow) 
bound  in  the  active  site  of  HuPONI.  The  proposed  catalytic  residues 
are  also  shown  in  color  (His  dyad  in  blue,  D269  in  red). 

catalytic  His  dyad,  implying  that  ester-  and  lactone-sub¬ 
strate  binding  adopt  a  similar  catalytic  mechanism  medi¬ 
ated  by  the  His  dyad  and  a  water  molecule. H>17  jn  con_ 
trast,  the  main  part  of  the  paraoxon  molecule  predomi¬ 
nantly  occupies  the  polar  binding  region  in  the  active 
site,  with  the  phosphate  group  pointing  towards  the  cata¬ 
lytic  residue  D269.  Interestingly,  MD  simulations  show 
that  the  distances  between  the  substrate  ester/lactone  car¬ 
bonyl  carbon  atom  and  HI  15  (s-N  of  imidazole)  remain 
stable  (~4  A),  whereas  the  distance  between  the  D269 
carboxyl  oxygen  atom  and  the  phosphorous  atom  of  par¬ 
aoxon  varied  between  4.8  and  3.5  A  (Supporting  Fig.  S4). 
In  Figure  10  the  predicted  binding  free  energy  decompo¬ 
sition  indicates  that  the  side  chain  of  D269  negatively 
impact  paraoxon  binding.  It  is  postulated  that  the  car¬ 
bonyl  group  of  D269  may  act  in  a  direct  nucleophilic 
attack  on  the  phosphorus  atom  of  paraoxon  in  the  cata¬ 


lytic  process  generating  a  high-energy  phosphor-enzyme 
intermediate.  The  unfavorable  interactions  of  the  D269 
side  chain  with  paraoxon  seem  to  be  compatible  with  a 
high-energy  intermediate  required  for  this  nucleophilic 
attack.  Probing  such  catalytic  mechanisms  involving 
bond-breaking  and  bond-formation,  however,  is  a  more 
challenging  problem  requiring  high  level  quantum  me¬ 
chanical  calculations. 

CONCLUSIONS 

Structural  information  describing  HuPONI  substrate 
interactions  is  of  primary  importance  in  the  development 
of  effective  HuPONI  variants  as  bioscavengers.  In  the 
present  study,  we  have  constructed  a  three-dimensional 
structural  model  of  HuPONI  and  probed  binding  inter¬ 
actions  of  HuPONI  with  various  substrates  (esters,  lac¬ 
tones,  and  phosphotriesters)  using  molecular  docking, 
MD  simulations,  and  binding  free  energy  calculations. 
We  have  delineated  dynamic  characteristics  of  the  active 
site  of  HuPONI  and  predicted  optimal  binding  confor¬ 
mations  associated  with  each  type  of  substrate.  The  pro¬ 
posed  binding  modes  and  a  number  of  key  residues  iden¬ 
tified  from  binding  free  energy  decomposition  supported 
the  conclusion  that  different  amino  acid  residues  in 
HuPONI  contribute  to  substrate  binding,  specificity,  and 
postulated  catalytic  mechanisms,  depending  on  the  class 
of  substrate  modeled.  Of  importance,  residue  Y71  is 
highly  associated  with  binding  of  all  three  modeled  sub¬ 
strates.  It  is  likely  that  this  residue  is  involved  in  the  cata¬ 
lytic  process,  serving  as  a  gate  in  the  active  site  to  facili¬ 
tate  substrate  recognition.  In  conjunction  with  these 
computational  studies,  we  are  currently  undertaking  site- 
directed  mutagenesis  studies  of  this  and  many  other 
amino  acid  residues.  Taken  together,  these  studies  shed 
light  on  the  structural  basis  for  substrate  interaction  and 
catalytic  mechanisms  associated  with  HuPONI,  and  pro¬ 
vide  a  rationale  for  further  engineering  of  this  enzyme  as 
a  catalytic  bioscavenger. 
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